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The medium modification of dihadron fragmentation functions in the deeply inelastic scattering 
(DIS) off a large nucleus is studied within the framework of the higher- twist expansion in the collinear 
factorization formalism. It is demonstrated that the modification due to multiple parton scattering 
in the nuclear medium is similar to that of the single hadron fragmentation function. However, the 
conditional distribution of the associated hadron given by the ratio of dihadron to single hadron 
fragmentation function shows only slight modification. The final results depend modestly on the 
nuclear density distribution. Comparisons with the experimental results on two hadron correlations 
as obtained by the HERMES collaboration at DESY are presented. 

PACS numbers: 12.38.Mh, ll.10.Wx, 25.75.Dw 

I. INTRODUCTION 

The modification of the properties of jets as they pass through dense matter has emerged as a new diagnostic tool 
in the study of the partonic structure of such an environment [l| . Such modification goes beyond a mere suppression 
of the leading hadron's multiplicity @, S, 0, H i, and could in principle be extended to include many particle 
observables [1, S E3, El, such as dihadron fragmentation functions. Defined as the two- hadron expectation values 
of partonic operators, these functions may be factorized [HI, [ld[ from the hard collisions and their evolution with the 
energy scale may be systematically studied in perturbative Quantum Chromodynamics (pQCD) [TEl, [l6l| . As the single 
hadron fragmentation, dihadron fragmentation functions in vacuum can be defined independently of the jet production 
processes. Once measured experimentally as semi-inclusive distributions of di-hadrons from one hard process, they 
can be used to predict dihadron distributions in other hard processes. The medium modification of the dihadron 
fragmentation functions in the multiple parton scattering picture will be dictated by jet-medium interaction as in 
the single hadron modification and therefore will provide an independent constraint on the jet transport parameters 
which in turn can provide important information on the properties of the dense matter. 

In both high-energy heavy-ion collisions and DIS off nuclei, two-hadron correlations have been measured (l7l. [Tsl. Il9| 
in addition to the single hadron spectra [2(], Hi], [HJ . Two-hadron correlations are measured as the conditional prob- 
ability distributions of a secondary hadron associated with a triggered hadron with given momentum. Therefore, 
they are given by the ratio of dihadron to single hadron fragmentation functions, while the single inclusive spectra 
are related to single hadron fragmentation function. The study of the medium modification of single and dihadron 
fragmentation functions in heavy-ion collisions involves a somewhat different kinematics than that in the DIS experi- 
ments and will be discussed separately. In this paper, the theoretical descripton of the dihadron correlations observed 
in DIS experiments off large nuclei will be presented. 

In conformity with the conventional picture of DIS on nucleons, the general process may be sketched as in Fig. [2] 
One of the quarks in the incoming nucleus receives a hard momentum transfer ~ Q 2 from the virtual photon. The 
struck quark then travels the remaning length of the nucleus while undergoing multiple scattering off soft gluons 
from various nucleons in its path. It escapes the nuclear environment and after some time, characterized by the 
confinement scale and dialated by its boost, undergoes the fragmentation process into a shower of hadrons. The 
presence of the hard scale (Q 2 ), provided by the momentum transfer between the lepton and the nucleus, allows for 
a factorized approach to the process of DIS off a nucleus, where the hard partonic parts can be calculated reliably 
within perturbative Quantum Chromodynamics (pQCD). At leading order, the differential hadronic cross section is 
estimated as a convolution of the quark distribution function in a nucleus with a leading order partonic cross section 
and the requisite fragmentation functions. The presence of the medium will be felt by the propagating quark through 
multiple final state scattering at next-to-leading twist level and the effects are in general suppressed by the hard scale 
Q 2 , While the final state interaction may be ignored in DIS off a nucleon, they play an important role in the DIS 
off a large nucleus where they are enhanced by the nuclear size [y^ 1 / 3 )]. This is due to the fact that the subsequent 
scattering with the gluon field (confined in nucleons) may occur at any of the nucleons which lie along the path of the 
struck quark. Such medium effects can be factoriscd from the leading order hard cross section and manifest themselves 
as medium modified structure [23| and fragmentation functions [|[ with the leading order (LO) hard partonic part 
remains unchanged. 
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In this paper, we will focus our attention on the medium modification of the fragmentation functions. The nuclear 
modification of the quark distribution functions will be included paramctrically. We will work within the formalism 
of higher-twist expansion in Ref. [23| to calculate the nuclear enhanced power corrections to the semi-inclusive cross 
section of DIS off a large nucleus. The modification of the single fragmentation functions within this formalism has 
been studied in Ref. [4J. In Sec. 2 we give a brief review of the formalism and isolate the contibutions generic to 
the modification of the dihadron fragmentation functions. The equations governing the medium modification of the 
dihadron fragmentation functions will be demonstrated to assume a form similar to that of the Dokshitzer-Gribov- 
Lipatov-Altarelli-Parisi (DGLAP) evolution equations 0, HH, [2(| of the dihadron fragmentation functions in the 
vacuum [ll|, This is quite similar to the case of the modification of the single fragmentation functions. The 
additional term as compared to the DGLAP evolution equation will be a twist-four matrix element indentical to the 
case of the single fragmentation functions. This matrix element will have to be evaluated in a nucleus of size A, 
and depends on the nuclear density distribution. In Sec. 3 the general methodology of the evaluation of this matrix 
element is outlined and its dependence on the nuclear density distribution is discussed. Numerical results for the 
medium modification will be presented in Sec. 4 for two different nuclear density distributions and are compared with 
the experimental data from the HERMES experiment. We summarize and present our conclusions in Sec. 5. 

II. FORMALISM 

The focus of this article is restricted to the semi-inclusive process of DIS off a nucleus where hadrons are detected 
in the final state. In particular, we will consider the following process, 

e(Li) + A(p) — » e(L 2 ) + Mpi) + ^(pa) + (1) 

when at least two hadrons are detected in the final state, where Li and L 2 represent the momentum of the incoming 
and outgoing leptons. The incoming nucleus of atomic mass A is endowed with a momentum Ap. In the final state, 
two hadrons (/ii,/i2) with momenta pi,p 2 are detected. 

The kinematics is defined in a frame as sketched in Fig. [TJ In such a frame, the virtual photon 7* and the nucleus 
have momentum four vectors q, Pa given as, 

P A =A\p+, 0,0,0], 

where we continue with the notation of Ref. [ll| of denoting four vectors in bold face. In this frame, the Bjorken 
variable is defined as xb = Q 2 /2p + q~ ■ 




FIG. 1: The lorentz frame of the process where a nucleon in a large nucleus is struck by a hard space-like photon. 
The differential cross section of the semi-inclusive process with two detected hadrons may be expressed as 

E L2 E pl E P2 da'^ M _ a 2 EM E Pl E p dW^ 
d 3 L 2 d 3 Pl d 3 p 2 2tts 4Q 4 ^ ' d?p x d?p 2 

where s — (p + Li) 2 is the total invariant mass of the lepton nucleon system and the leptonic tensor is, 



Li 



-Q 

2q- 



q-,0,0 



(3) 
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In the notation used in this paper, \A,p) represents the initial state of an incoming nucleus with A nucleons each 
carrying momentum p. The hadronic part of the final state with two detected hadrons is represented as \X,hi,h 2 ). 
As a result The semi-inclusive hadronic tensor may be defined as 

Epi Ep2 d?p x d?p 2 = 2n5i - q + P A-PX-P1-P2) 

x (A,p\J^(0)\X,h 1 ,h 2 )(X,h 1 ,h 2 \r(Q)\A,p), (4) 

where the sum QZx) runs over au possible hadronic states and J M is the hadronic electromagnetic current (J' 1 = 
eqtpq'y^i^q). It is understood that the quark wavefunctions are written in the Heisenberg picture. The leptonic tensor 
will not be discussed further. The focus in the remaining of this paper will be exclusively on the hadronic tensor. 
This tensor will be expanded order by order at the leading log. The leading and ncxt-to-leading twist contributions 
that are nuclear enhanced will be isolated. 



A. Leading twist 

We start by noting that within the kinematics chosen there are two sets of final states in the opposite directions: 
one set \S\) is in the direction of the nuclear beam and consists mostly of its remnants, while the other \S 2 ) is in the 
direction of the hard photon and consists of a jet of hadrons. At leading twist, there is minimal overlap between these 
two sets and this leads to a simplified expression for the integrated hadronic tensor, 



d 3 p! d 3 p 2 
(2n) 3 2E 1 (2n) 3 2E 2 



J2 (27r) 4 <5 4 (q + P A - PSl - Ps 2 Pi P2) 



X p h\Jl<2 S1S2 

x (Mp|^(0)|5 1 )7' 1 (0|^(0)|5 2 p 1 p 2 }(^ 1 p 2 |^(0)|0}7' y (5 1 |^(0)|pAp). 



(5) 

In the above equation, we include an average over the intial spin states X. p of the incoming nucleus, in the following 
this will be implied. We have also simplified \A;p) as simply \p); the presense of A should always be understood. 
The above equation is restricted to the case of a quark inside the target undergoing a hard collision with the virtual 
photon. The case where an antiquark is struck is analogous and will not be discussed. The above is the leading 
order term in the usual expansion of the quark wavefunction operators in the interaction picture. In the rest of this 
section, the factorization of the dihadron fragmentation function from the hard part and the structure function will 
be discussed. It will be demonstrated that no new assumptions need to be invoked. The short hand notation 



11 h-2 



E 

hi ,fi2 



pi 



P2 



(2tt) 3 2E 1 {2it) 3 E 2 



will be used henceforth. As in Ref. we use the notation that p\ + p 2 = ph- 

The ^-function can be espressed as an integral over four space-time, where phase factors involving the proton's 
momentum and that of its remnants may be absorbed into the matrix elements of the quark wavefunction operator 
on the proton state. Under a reorganization of the spinor matrix elements and introducing the momentum fractions 
£, one obtains the hadronic tensor as 



W<™ = ^ J d 4 xe ix(q - ps 2-P^ J d 4 kd(6 4 (p S2 +Ph-k-q)5(( 

^(5 1 K<o)b)(piv;(x)^ 1 )7 Ai <o^(o)i^ 1 p 2 )(5 2 p 1 p 2 |v;(o)|o)7 !/ 



p+) 



x Tr 



(G) 



In the above, the four-vector k has been introduced which represents the momentum of the struck quark inside the 
nucleus. Using the four (5-function, one may re-express exp[ix • (q — ps 2 — p^)] as simply exp[— ix • k]. Taking the 
Fourier transform of the <5-function over fc, leads to the expression, 



x Tr 



' V^- j;x ' k <5(C - ^) 



(SMmp)W(x)\Si)r / d 4 ye 



-iy-(Cp+q) 



(0|^(0)|5 2 p 1 p 2 )(5 2 p 1 p 2 |^(0)|0)7 l/ 



(7) 
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Invoking the collinear approximation that the (+)-components of the target momenta p + , k + are much larger than 
their (_L)-componcnts and (— )-components, one may approximate the leading twist contribution of the x dependent 
matrix element as 



f(C ' p) = J ^d*xe-™- k 6((-^)Y,(Si\m\p)(pmx)\Si)^p(0- (8) 

Si 

Where, /(£) is identified as the unpolarized quark distribution function in a nucleus, which now has the obvious 
expression, 

/(c) = J ^^xe-^Hic-^ipmx^mip)- (9) 

Thus, the initial quark distribution function may be factorized from the hard cross section and the final fragmenta- 
tion into hadrons under the collinear approiximation. Factorization of the dihadron fragmentation function from the 
hard part requires the use of the collinear approximation in the opposite direction of the produced hadronic jet. One 
introduces the quark momentum in the direction of the hadronic jet as 



J d 4 /<5 4 (l-(Cp + q)), 



as well as the momentum fractions of the two detected hadrons as 

!i — -rr) 5 {%2 — tt 



We continue with the notation of Ref. 11| where z = z\ + Z2- As a result, the simplified expression for the differential 
hadronic tensor for the production of two hadrons with momentum fractions z\, Z2 in the direction of the quark jet is 
obtained as 




FIG. 2: The Lowest order contribution to . 

The double integral over the momenta of the detected hadrons may be re-cxprcsscd as an integral over the sum of 
the two momenta and the difference. The functions S(p~^z — q~ — (p~)6 2 (ph±_) may be invoked to restrict the sum of 
the momenta of the detected hadrons. While in the limit afforded by the collinear approximation, the (^-function over 
the (+)-componncts of the momenta yields, 



'l^-.*-tf*j=«(^=-<» 
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In the above, x B is the Bjorken scalin g va riable. The two 5-functions over the two momentum fractions z\, z 2 may 
be readjusted as demonstrated in Ref. [ll[ which leads to the final factorized expression for the differential hadronic 
tensor as, 



= / dx B f(x B )H^D^ h *(z u z 2 ), 



dz\dz2 

where, H^ flL ' is the hard part of the quark-photon scattering i.e. 



(11) 



and Dq 1,h2 (zi, z?) is the dihadron fragmentation function, 



(12) 



D 



hi } h$ 



d 2 q± 



(2tt) 4 { I 



4ziz 2 J 4(2tt) 3 
x (0\Mx)\ Pl , P2 ,S - 2){p uP2 , S - 2|^(0)|0) 
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S-2 



(13) 



In the above expression, q± = pu_ — p2j_. This expression is identical to that presented in Ref. [ll[ with the obvious 
switch in the direction of the final momentum from the (+) direction to the (— ) direction. Also, Eq. (|12p may be 
directly generalized from the expression for single inclusive hadron production presented in Ref. The resulting 
picture in terms of Feynman diagrams is thus that of Fig. Hence, it may be stated that the factorization of the 
dihadron fragmentation from the LO diagram in DIS may be easily generalized from the factorization of the dihadron 
fragmentation function in e + e~ annihilation. 

The method of factorization of the dihadron fragmentation function followed here may be easily repeated in the 
case at next-to-leading order to which we now generalize. The case at higher twist will be discussed in the next 
subsection. At next-leading-order in the strong coupling constant a s , the dominant contribution at the leading log 
level (in n ■ A = gauge) comes from final state gluon radiation as shown in Fig. [3J Specifically, these large logarithms 
result as the transverse momentum of the radiated gluon vanishes (i.e., lj_ — > 0). In this limit, the two detected 
hadrons may materialize in the fragmentation of the struck quark or from the radiated gluon. It is also possible that 
one hadron originates in the fragmentation of the quark while the other originates in the fragmentation of the gluon. 
In all these cases, the fragmentation functions may be factorised as a convolution with the hard cross section and 
a set of splitting functions, which encode the probability of a quark radiating a gluon which carries away a certain 
fraction of its forward momentum. 

No doubt, these contributions result in the Dokshitzer-Gribov-Lipatov-Altarelli-Parisi (DGLAP) evolution of the 
dihadron fragmentation functions. These have been discussed at length in Refs. [ll], E3 f° r the case of jets in 
e + e~collisions. Using the methodology in these references and the preceding discussion at leading order, the leading 
log contributions to the hadronic tensor, as the radiated gluon becomes collinear to the struck quark, of the semi- 
inclusive DIS process may be expressed as 




In the above equation, Hj^J represents the leading order hard cross section in Eq. (TT2"|) . The quantities P q ^ qg (z) 
represent the probabilities that the struck quark radiated a gluon and retained a fraction z of its energy. They have 
the well known expressions 0, HH, Ull , 
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Pq ^ qg ( z ) = C F (^Pj . (15) 

The subscript '+' indicates that the negative contribution from a virtual correction has been added within the 
splitting function. The splitting function in the second line of Eq. (fT4)l . P q ^ qg , is identical to the above equation 
except that it lacks the negative virtual correction. The other splitting function P q ^ gq (z) is obtained simply as 
P g _, 9s (l — z) and also does not admit a virtual correction. 

The gluon dihadron fragmentation function D g lh2 (zi, Z2) is obtained as 



sab Jfiv n\ 

x Y / (0\A;(x)\ Pu p 2 ,S-2)(p 1 ,p 2 ,S-2\Al(0)\0)°-—^l, (16) 

S-2 

where d^il) is the gluon's polarization tensor in the light-cone gauge and sum over the color indices of the gluon 
field is implied. The meaning of the various momenta and momentum fractions in Eq. (|16[) is identical to the case 
of the quark fragmentation function. The single fragmentation functions D q (z), D g (z) have the usual expressions 

The NLO modification at leading twist presented above suggests the formulation of an effective next-leading-order 
and leading twist hard part. This includes the imaginary part of the photon nucleus forward scattering amplitude at 
order a s . In particular, the hard part represents the process of a quark in the nucleus being struck by a hard photon 
into an off-shell intermediate state which then radiates a gluon, 



H^( x ,p, q ,z) = H^( x ,p, q )£ <^7 F i±i (17) 

This hard part is represented by the diagram in Fig. [3l One can in principle deal solely with this hard part. The 
requisite double and single fragmentation functions may be convoluted with this hard part as a separate final step. 
This allows one to focus solely on the derivation of the hard part at the leading twist. Note that this hard part 
is identical to that derived for the case of the modification of the single fragmentation functions in Ref. 0]. There 
remains the question as to whether such a separation may be achieved at higher twist. This is the subject of the next 
subsection. 




FIG. 3: The hard part of the next-to-leading order contribution to W 



B. Higher Twist 



In the calculation of cross sections which involve a hard scale, it often suffices to simply calculate the leading 
twist part. Higher twist contributions are suppressed by powers of the hard scale. It was demonstrated in Ref. [23j |. 
that a class of higher twist operators are enhanced by the size of the target, A 1 / 3 , where A is the atomic mass. 
Specifically these contributions arise from the multiple scattering encountered by the struck quark off the soft gluons 
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as it bore through the nucleus. Such contributions are also important in the medium modification of the single hadron 
fragmentation functions In this section we isolate the hard part of the diagrams involved in the modification of 
the dihadron fragmentation functions. 

In the previous subsection, we demonstrated the factorization of the fragmentation and parton distribution functions 
from the hard scattering cross section at leading and next-to- leading order and at leading twist. It was demonstrated 
that the hard part is the same as in the case of the single fragmentation functions. The focus will now be on a similar 
proof at next-to- leading twist. At leading order the diagrams to be considered are illustrated those in Fig. [U 




FIG. 4: The leading order and next to leading twist contribution to W^" ' . 



The reader will immediately note that the procedure involved in the factorization of the dihadron fragmentation 
function from the hard part is entirely similar to that at leading twist, and also to the case of the single fragmentation 
functions. The resulting hard part has no radiated gluon and therefore no medium modification of the fragmentation 
functions. The sole effect of the double scattering indicated by the diagrams in Fig. [3] is the appearence of a net 
transverse momentum k± in the out going jet brought in by the average transverse momentum of the gluons. In the 
limit of vanishing k± , the contribution of the above diagrams is an overall Eikonal phase that will be absorded into a 
gauge invariant definintion of the quark distribution function. 

Medium modifications to the semi-inclusive cross section at higher twist in the limit of very small transverse 
momentum of the soft gluons in a nucleus emanate specifically at next to leading order. The contributions are 
illustrated by diagrams in Fig. [5J Both diagrams represent the case where the struck quark undergoes scattering off a 
soft gluon in the nucleus and also radiates a fraction of its forward momentum into a collinear gluon. Both the quark 
and gluon then exit the medium and fragment. The diagram on the left corresponds to the case where both detected 
hadrons emanate from the fragmentation of the quark, while the diagram on the right corresponds to the case where 
one of the detected hadrons originates in the fragmentation of the quark and the other in the fragmentation of the 
gluon. In Fig. Owe two diagrams that arise at next-to-leading order and at next-to-leading twist in the evaluation of 
the medium modification of the dihadron fragmentation function. The dominant contributions, however, comes from 
gluon-gluon secondary scattering 0]. 




FIG. 5: A contribution at next-to-leading order and next-to- leading twist to W^ v . 
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As an example, we focus on the diagram on the right: the case of two independent fragmentations. The factorization 
of the dihadron fragmentation function from a quark line, as is the case for the diagram on the left in Fig. [5] is a 
simpler generalization from the case of the single inclusive cross section. Thus, the case of two single fragmentations 
represents the new addition to the type of diagrams involved in the calculation of the medium modification of the 
dihadron fragmentation functions. It will be demonstrated that the next-to-leading twist hard part, which arises 
from the evaluation of this contribution, is no different than the hard part, at next-to-leading twist, in the case of 
the single fragmentation functions. The difference lies in the convolutions involved: besides the usual convolution 
with the quark and gluon dihadron fragmentation functions, one also needs to convolute with a product of two single 
fragmentation functions as in Eq. (fT4|) . 

As an illustraion, the contribution to the hadronic tensor from the right-hand diagram in Fig. O may written in 
general as 



* <*>Ht<0)|0> J ||^^ 7 ^ e -.,-<- + .>e-<-<'.-^>r(0W(O)| S2 > 1 ) 

x (sjpiWojloh"/ ^«' yrll ' _k, F f7;^^7p^7»(Sil^(»)|Jf)i;<x|*(0)|p)- (is) 

In the above equation, the gluon momentum / = P S b + p 2 and quark momentum l q = Ps» + pi arc implicit. The sums 
Sx'Sy represent intermediate states assumed by the target after the initial quark is struck by the hard photon. 
The final state sums are 5Zg in the direction of the target, X^s 2 m the direction of the photon. The latter is split 
into two parts S 2 emanating from the fragmentation of the quark and 5*2 emanting from the fragmentation of the 
gluon. The momenta k, k 2 represent the momentum of the quark between the radiation of the gluon and scattering 
off the soft gluon field. 

We introduce the following factor of unity: 

/ d ^ / J J ^J^-^ J dz[5 (4 _ |) / dz > S (4 - fj , 

followed by a reorganization of the spinor matrix structure as in Eq. ([6]) which leads to the trace over spinor indices. 
The introduction of the various integrals and 5-functions allow for a simple procedure of indentifying and extracting 
the matrix elements in the fragmentation functions. Using the (^-functions over the momentum fractions, one may 
re-express the partonic momenta l q , I in terms of the rescaled hadronic momenta p\j 'z[,p 2 / ' z' 2 . As a result the quark 
fragmctation function may be isolated as, 



. i 3 



dH. 



2 J (2tt) 

and the gluon fragmentation function as, 



2pi 



-W(a:,)|S2pi><S2pi|V(0)|0) 



(19) 



D h g 2 (z' 2 ) = \\ (|^(^-p) / ^e^^OlA^x^l^XSSpal^^lO)^^. (20) 
Utilizing the above relations, one obtains the differential hadronic tensor as 



dW"'" = dz[dz' 2 e 2 g 4 Tr[t d t c t a t b ] / d 4 yd i y l d i y 2 



J tll,h 2 L g 



d i kd i k 2 fa + / 
(2tt) 8 {k 2 + l) 2 -ie 



r y> ■ 1 ■■■■ ■. ■ 

' k 2 -ie' 2 ' k 2 +ie (k + I) 2 + ie 



d af> 5 



iy<q-k 2 -^-) <»-(fa-2f) iyf(^--k) DJz[) D g {z' 2 ) 
x e " 2 e 'I e -i q y . 

(zir (4) 



(21) 
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The factorization of the fragmentation functions from the hard part is completed by reexpressing the hadronic 
momentum integrals in terms of their rescaled partonic counterparts i.e., d 3 pi ~ d 3 l q (z[) 3 and d 3 p 2 — d 3 l(z' 2 ) 3 . We 
also introduce the variable z, the fraction of momentum that remains in the quark after it has radiated a gluon. In 
terms of z the hadronic momenta are rescaled to z\ = z[z and z 2 = z' 2 (l — z). Stipulating that the matrix elements 
of the initial states do not depend strongly on y^, y\^ + y 2 ^, the integrals over the transverse coordinates may be 
completed yielding (5-functions over the transverse momenta. The matrix elements involving the initial states are 
dominated by the (-) component of their coordinates. Integrating over these transverse momenta under the influence 
of the above mentioned 5-functions, the double differential hadronic tensor may be expressed as 

x ie»'*'(H^(»-)T + -«' t W,nM" , '(»i'.0)*<0)|P). (22) 

The reader will note that the double differential cross section has been expressed as a convolution of an initial 
state quark-gluon correlation in the incoming nucleon(nucleus), a hard scattering piece (H^(p, q, z)) and final state 
quark and gluon fragmentation functions (D a (x), D a (x)). Isolation of the initial state piece involves the well known 
decomposition of the gluon vector potentials |23| , 



where, 



A% ~av, CT A Q</ +p CT — -, (23) 
p ■ n 



P a np + ppn a 

J a,/3 — 9aj3 — 



p ■ n 



Within the given kinematics, the contribution from the second term in the Eq. (|23j) far out weighs that from the first 
term. This is followed by the usual approximation of isolating the leading twist piece of the quark gluon correlation 
function. The hard part in Eq. (f2"2"| is obtained as 



„ m , . ,Tr\t d t a t a t' lr 
H<U(p,q,z) = g 4 L , Tr 



fr_^ v I dx 1 dx 2 dxidk\ dH $+^2^+^ A + %a t> 



(2tt) 5 (2tt) 4 



x 7 



x 7" 



0. 



2q ■ pz(xs - x± + x L + ie) 

A + i> 



i>{A + (x-l + x 2 ) jp— /+ fc±) i) 



2q ■ p(xb - X4 + ie) 

A + x i i>- / 



2q ■ pz(xb — X\ + x L - ie) 



2q -p(x B - xi - ie) 



d a pS 1 — z — 



(24) 



In the above equation, x±,x 2 ,X4 represent the forward momentum fractions of the incoming quark and gluon and 
outgoing quark in Fig. O The reader will also note that the hard part presented above in the evaluation of the double 
inclusive cross section is identical to the hard part in the evaluation of the single inclusive cross section Q ■ 

The remaining analysis of the cross section involves expanding the hard part as a Taylor expansion around the point 
k±_ = 0. The leading term in such an expansion yields the first term in the Eikonal expansion of the intermediate 
quark propagator in Fig. [3l The second term which involves a linear derivative of the transverse momentum is zero 
identially for unpolarizcd targets. The third term survives and represents the first contibution to the nuclear enhanced 
higher twist contributions to semi-inclusive deep inelastic scattering. The evaluation of the higher twist contribution, 
factorized from the fragmentation functions, is in every way identical to that for the single inclusive cross section Q. 
Extracting the leading order hard part H 1 -^ from and reorganizing the double inclusive cross section in the form 
of Eq. dHJ), i.e., 



~dz~^z~ 2 = J dx Bf(x B )H^D h q ^ h2 (z 1 ,z 2 ), 
one obtains the medium modified dihadron fragmentation function Dg u "(zi, z 2 ). 



(25) 
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The modification of the fragmentation function in this sense depends on the scattering of the struck quark off the 
gluons in the medium and is thus dependent on the initial state. Following the methods employed in the computation of 
the medium modification of the single fragmentation functions, one obtains the medium modification of the dihadron 
fragmentation functions as the sum of higher twist contributions extracted from the sum of multiple diagrams at 
next-to-leading order. The presence of multiple fragmentation options allows for a greater number of contibutions to 
the equations governing the medium modification of the dihadron fragmentation functions. These may however be 
grouped together in a similar fashion to the vacuum evolution contributions to the dihadron fragmentation functions. 
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2 
_L 

I 2 2tt 



dl\ a 



dz 



h *( Zl ,z 2 ) = D h q ^( Zl ,z 2 )- 

+ AP q ^ gq (z,x B ,x L ,ll)D^ I' 

\ z z / > J Z1 



AP q ^ qg (z,x B ,x L ,li)D^ 
l ~*> dz 



h\ ,h 2 



Zl £2 

Z ' y 



z(l 



AP q ^ qg (z,XB,XL,l 2 L ) 



x D 



■t) D 



Z-2 



1 



{hi 



(26) 



qg , AP q ^ gq represent the 



In the above z^ = z\ + z 2 , and the switch (hi — > h 2 ) is only meant for the last term. AP q 
medium modified splitting functions where a momentum fraction z is left in the quark and the gluon respectively. Their 
expressions are identical to the modified splitting functions derived in Ref. In the above, xl = l±/(^P + q~ z(l —z)), 
where the radiated gluon or quark carries away a transverse momentum 

As in the case for the evolution in the vacuum (Eq. (TH)) ). the splitting function P q ^ qg 
counterpart and is given as 



= 9\t[P q ^ qg ] has no virtual 



l + z 2 C A 2^aT^ g (x,x L ) 

i-z {ii + (ki))N c mx, l p I r 



^■Pq^qg — -, _ /,2 , /7„2\\Ar tA(„ ,.2V (27) 



In the above equation, Ca = 3, N c = 3. The scale \i\ represents the scale at which the quark distribution functions are 
factorized from the hard cross section at leading order and leading twist. The mean transverse momentum of the soft 
gluons is represented by the factor (fc^). The term represents the quark gluon correlation in the nuclear medium 
and includes contributions from multiple higher twist diagrams. The formal derivation of T qg , identical to the case 
for the single fragmentation functions is outlined in Ref. [J] . The final result which includes squares of amplitudes of 
soft-hard scatterings, hard double scattering and their interferences results in the simplified form, 

T*(x B , x L ) = J ^dy^dy 2 e^ XB+XL ^ +v '(l-e-' lXLP+y 2)(l~e-~ lXLP+ ^- y ^) 

x \ip\mi + F+{y 2 )F +a {y^(y-)\p)6{-y-)e{y 2 -y^). (28) 

Up to this point the state \p) was generically referred to as that of a nucleus with momentum Ap. Indeed, the isolation 
of the higher twist piece in DIS so far has been independent of the content of this state. If \p) represented a nucleon 
with momentum p then the above equation would represent the medium modification factor for DIS off a nucleon. The 
enhancement of such higher twist objects in extended nuclear media will be the subject of the subsequent sections. 
This will entail the decomposition of the nuclear medium into nuclcons and will directly involve the nuclear density 
distribution. The evaluation of T qg will expose the eventual nuclear size enhancement and as a result justify the 
incorporation of the higher twist contributions in the evaluation of inclusive hadron production in the DIS off a large 
nucleus. 



III. EVALUATION OF THE NUCLEAR MODIFICATION FACTOR 



In the preceding section, the higher-twist contribution to the inclusive hadron production in deep-inelastic scattering 
off a nucleus was expressed as the convolution of the quark distribution function in a nucleus, a hard LO photon-quark 
scattering cross section and a nuclear modified fragmentation function. The equations governing the modification of 
the dihadron fragmentation function in medium (Eq. (|2"E|) ) depend on the nuclear modification factor (Eq. (l28l) ). 
This quantity essentially represents a quark-gluon correlation function within the nucleus. The evaluation of this 
correlation function is the subject of this section. 
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The essential kinematics for the process of DIS off a large nucleus has been outlined in Sec. II. A virtual photon 
with momentum q strikes a large nucleus with atomic mass A with a momentum Ap. In the evaluation of the 
multiple scattering of the struck quark in a large nucleus we will invoke the convolution model, which decomposes 
the deep-inelastic scattering off a nucleus in terms of the scattering off its constitutient nucleons. In this paper, we 
will follow the version of the convolution model outlined in Ref. [29j . The evaluation of the quark-gluon correlation 
function which arises as the struck quark traverses such a nucleus commences with the decomposition of the nuclear 
ket \p; A) in terms of nucleons, 

M> =[f[ S#T*(fe})lfe})(2-) 3 2 P + 5 3 (^ ? ? J - Ap). (29) 

i—l v / i 

In the above equation, $({pi}) represents the nuclear wavefunction; the ket = \pi,p 2 , . . . ,p A ) represents a 

particular ket where the nucleons 1 . . . A assume momenta p\, . . . ,p A - The overall three momentum conservation is 
enforced via the (5-function. In the notation employed in this article, a single nucleon state is normalized as 

(p*bj> = 2p+(27r) 3 5 a (p l - ft), (30) 
as a result the n-nucleon state is normalized as 

n 

({Pi}\{Pj}) =Y[2pt(2nf6 3 (p i -p j ). (31) 

i=i 

Defining the nuclear momentum Pa = Ap, the normalization of the nuclear state may be stipulated as, 

{p- A\p'; A) = 2p+(2n) 3 S 3 (P A - P' A ). (32) 

Substituting the expression for \p; A) from Eq. (|29[) into Eq. (|32[) leads to the normalization condition on the nuclear 
wavefunction (where we use the simplified notation d]jf = ]jf =1 d 3 p l 9{p+) /[(27r) 3 2p+]): 

A 

/ d^ttp^Mip^n^p+S 3 ^^ - Pa) = 1. (33) 

^ 1 i 

Halting one of the momentum integrations on the l.h.s. of the above equation leads to the momentum space 
one-nucleon density, which represents the probability of finding a nucleon with momentum k, 

. A A A 

p(k)= dl[\<P(k,{p l })\ 2 (2n) 3 2p + 5(Ap + -k+ -J2p t )s 2 {k±-J2P^)' ( 34 ) 

■* 2 i=2 2 

where, wc have decomposed the 3-5 function into its longitudinal and transverse components. Halting two integrations 
leads to the two-nucleon momentum correlator fe, A), which contains information regarding the sharing of 
momentum between nucleons in a nucleus, 

A 

p(k!,h,A) = d]J^*(ki + A/2, k 2 - A/2, {pi})^(k!- A/2, k 2 + A/2, { Pi }) 

3 

A A 

x (2TT) 3 2p + d(Ap + -k+ -k+ -Y / P l )s 2 (k h± + k 2 , ± -Y / &±)- (35) 



i—3 i—3 

Both these and the higher correlations obey the obvious sum rules, e.g., 



J dkp(k) = 1, J dkidk 2 p{ki,k 2 ,0) = 1, (36) 
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where we have used the short hand dk = d 3 k9(k + )/[(2ir) 3 2k + ]. 

The matrix element included in the definition of the modification factor in Eq. (|28[) involves interferences between 
the different processes which contribute to different phase factors in Eq. (|28|) . The interference between the various 
processes is thus dependent on these factors and ultimately on the locations y,y\,yi where the various scatterings 
have taken place. Hence, the nuclear wavefunctions discussed above need to be transformed to position space. In the 
notation adopted in this paper, the position bases demonstrate the following outer product and normalization, 



n^|{^}>({z 4 }| = 1, ( Xi \ Pi ) = e&^yfitf. (37) 

i 

Within this basis, the nuclear ket may be expressed as \A) = J Yli d 3 Xi$(xi)\{xi}). To simplify the expessions, we 
use the over all (5-function in Eq. (|29[) to carry out the integration of the A th momentum. There are now simply A — 1 
nucleon momenta to be dertcrmincd to isolate as kiven ket, the wavefunction provides the probability amplitude for 
a given choice of A — 1 momenta such that one of them is specified. The momentum of the final ket pa is such that 
the momentum of the nucleus remains Pa i.e., p\ = ^a ~ 2* 1 Pi • Introducing a complete basis of position states 
on the r.h.s. of Eq. (|29[) . one can reexpress the nuclear state in position space as 



.A-l „A-1 



\p;A) = fU d 3 x t [ [] dp. l ^dmp l },p A )e^P—2p + \{x l }). (38) 
J i=i JZPa 



By mere inspection, the nuclear wave function in position space may be expressed as 



*({*<})= fl[dp l ^^Mp l ,pA)e^t^2p + , (39) 



where, once again, it is understood that pa = Pa — Tlit = i Pi- As a result the wavefunction in position space has the 
right translational properties, i.e., 

$({£, + Ax}) = e lPAAx <S>{{x,}). 

The inverse transform is obtained as 



A-l 

2P+-2]Tp+ (40) 

i 

Armed with the decompositions above for the nuclear wavefunction, one can now evaluatr the nuclear modification 
factor in Eq. (|2"8f . The primary object to be evaluated in Eq. (|28p is the matrix element 

M = (p; A\^(0h + F, •+(y^)F +a (y^)^(y~)\p;A). 

The quark and gluon operators carry color. However, the expectation of this operator is evaluated in a nucleus which 
is color neutral. In the convolution model outlined above the nucleus is decomposed into its constituent color neutral 
nucleons. The various position variables 0, y~ , y^, 2/2" represent the locations in the nucleus where the quark was 
struck by the hard photon and where the soft gluon was encountered by the struck quark. As the quark and gluon 
carry color in different representations of SU(3), both quark operators must be confined within the same nucleon, 
as must both gluon operators. It is possible that all four operators are confined within the same nucleon. This is 
definitely the case for the higher twist contribution to DIS off a nucleon. In a nucleus, the two quarks operators and 
the two gluon operators may be confined to two different nucleons. In the limit in which the struck quark carries a 
large forward momentum x\p + and the soft gluon carries a very small transverse momentum k±, the struck quark 
follows an almost linear trajectory as it burrows through the nucleus. As a result the struck gluon may be located in 
any of the nucleons that lie directly behind the nucleon which is struck by the photon. This leads to an enhancement 
of the order of one nuclear dimension i.e., ~ A 1 / 3 . In this and the next section, this factor of A 1 / 3 will be explicitly 
extracted. This nuclear size enhancement, aposteriori, justifies the following approximation for the matrix clement, 



*({Pi» = — / ] 



A-l 



1 -. _ 

Pi 



*«*}) 



\ 
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* = (2) / n* 



Of\ /" d 3 k l d 3 k 2 d 3 A6{kf + A+ /2)(9(fc+ - A+ /2)8{k+ + A+ /2)6>(fc+ - A+ /2) 0(p+) 

^ a -1-2 / +2 a , 2 9„+ 



(2*)*2*{kf - A+ 2 /4)(fc+ z - A+ 2 /4) W A 
x $*(fex - A/2, k 2 + A/2, {pi}^ 1 )*^! + A/2, k 2 - A/2, { K }^- 1 )(2p+) 2 

x (h - A/2|^(0) T +V(y-)|fci + A/2)(fc 2 + A/2|F+(y 2 -)^ + (yr)|fc 2 - A/2). (41) 

In the above equation, we have simply reexpressed the nuclear ket in terms of nucleon kets in momentum space as 
required by Eq. (f2T))) . The set of states {Pi}^ -1 represents 3 to A — 1 nuclcons, the momenta of the first two nucleons 
are k\ ± A/2 and k 2 ± A/2. Due to color confinement the two quark operators in M act on the same nucleon, as do 
the two gluon operators. The factor of (^) originates in the number of ways the two nucleons may be chosen. The 
remaining A — 2 nuclcons are unaffected by the operators and lead to the appearence of A — 2 momentum (5-functions 
as in Eq. (|3U)) . Integrating over these momenta removes A — 2 of the 2 A — 2 integrals leaving the A integrals as in 
the equation above. The momenta of the two nucleons on the r.h.s. of the Eq. (|4"Tj) which form the bras and kets of 
the two partonic operators are such that the sum of the momenta of the two nucleons that form the bras are exactly 
equal to the sum of the momenta of the two nucleons that form the kets and each is equal to k\ + k 2 . This mismatch 
in momentum is the reason for the appearence of the two off-forward parton distributions [3(| which will result from 
the fourier transform of the following shifted matrix elements of the partonic operators, 

(h - A/2|^(0) 7 + ^(y-)|fc 1 + A/2)(fc 2 + A/2\F+(y-)F-+(y-)\k 2 - A/2). 



Utilizing the expressions afforded by the preceding discussion, the matrix element in Eq. (|41|) may be transformed 
to position space. The A — 2 6- functions which contain all but the mismatched momenta i.e. p^, . . .p\ are Fourier 
transformed, 



/ r j 7 . piZiPj 
I™ < 42 > 

The A — 2 momentum integrations can now be carried out setting the position of the \ th nucleon in $* to x\ = Xi + z%. 
Utilizing the fact that the nuclear wavefunctions are rapidly dropping functions as any one of the coordinates tends 
to infinity allows the evaluation of the A — 2 countour integrations over the variables z,, yielding the expression 



M 



A} f d 3 k 1 d 3 k 2 d 3 A6(kf + A+/2)0(k+ - A+/2)0(k 



A+/2)0(ft+ - A+/2) 



(27r) 9 2 2 y (fc+" - A+74)(fc+ - A+74) 
x / d 3 x 1 d 3 x 2 d 3 x' 1 d 3 x' 2 $*(x' 1 ,x' 2 ,{x l })<f>(xi,x 2 ,{x l }) 



- i(h + A/2) • Xl - i(k 2 - A/2) • x 2 + i(h - A/2) • x[ + i(k 2 - A/2) • x 2 



x exp 



x (h - A/2|^(0) 7 +V(jr)|fci + A/2)(fc 2 + A/2\F+{y-)F°+{y-)\k 2 - A/2). 



(43) 



^From Eq. (|28p . the complete nuclear modification factor is composed of four terms, each of which contains the 
matrix element expressed in Eq. (|43[) . The four terms represent the squares of amplitudes and the interferences 
between double hard and hard soft scattering. Each of these terms may be generically expressed as 



T = |^dyr<foae W "~^ flP+( ^^^^ ( 44 ) 

where, xa, xb, %c represent different combinations of momentum fractions, and M(y — , y^, y 2 ) represents the same 
matrix elements as those of Eqs. (|4"Tj) and (|43[) with the dependence on the partonic position variables shown explicitly. 
As the reader will note, the full form of T is essentially a convolution between two factors one that contains an integral 
over nucleon position variables and another that contains an integral over partonic position variables. These two terms 
are convoluted through the nucleon momenta k\, fc 2 , A i.e. 



T= d 3 k 1 d 3 k 2 d 3 A[N(k 1} k 2 , K)P(ki,k 2 , A)] , 



(45) 



14 



where, 



and 



N(k!,k 2 , A) = J d 3 xid 3 x 2 d 3 x[d 3 x 2 N(k 1 ,k 2 ,A,x l7 x 2 ,x' 1 ,x 2 ) (46) 

P(ki,k 2 ,A) = J dy~dy^dy 2 'P(k 1 ,k 2 ,A,y-,y^,y 2 -). (47) 

New partonic position variables with unit Jacobian vT = y-j" — and = (y^ + y 2 )/2 arc introduced. In terms 
of these the partonic matrix elements of the gluon operators may be simplified as 

(k 2 + A/2\F+(y 2 )F^(y^)\k 2 A/2) = (k 2 + A/2\F+ \k 2 - A/2)e^ +w ~ . (48) 

Similarly, the ^-functions over the partonic locations may be reexpressed as 

f dkdl p - ik ( w ~ -«~ / 2 ) p—U(v~— w~-u~/2) 

0{-y 2 )6{y- - yr) - '(«"/2 - - to" - «"/2) = / 7^^— 7" ■ (49) 

Integrating out w~ one obtains the constraint 27ri5(A + — (k + I — xcp + ))- This 5-function is used to perform the 
integral over the longitudinal off-set in the momenta i.e., A + . It has already been stated that the matrix elements of 
the partonic operators in Eqs. (I41I43[) are. under the action of Fourier transforms, related to the off-forward parton 
distributions (OFPD) |30| . However, as was pointed out in Ref. [3l|, MIT bag models of the OFPD's suggest that the 
variation with A is rather slow and thus we will assume that the OFPD's remain more or less constant with A. It will 
also be assumed that the matrix elements are almost independent of the transverse momentum of the nucleons i.e., 
ki±, k 2 j_. This allows the integration over k±±, k 2 ± and Aj_ which leads to three sets of (5-functions over the transverse 
components of the nucleons involved: S 2 (xi± — x' 1± ) , 5 2 (x 2 ± — x 2J _) and S 2 (xi± — x 2 ±). Thus all the transverse 
positions of the nucleons are set to be the same. As a result, the simple linear trajectory of the propagation of the 
struck quark as presented in Sees. II & III is reinstated. The representative matrix element, with the incorporation 
of the above simplifications, assumes the form 



T 



dkdl f dk+dk+ 9(k+ + A+ /2)6{k+ - A+ /2)6(k+ + A+ /2)9{k+ - A+ /2) 



^ 4 v /(fc 1 - 2 -A+ 2 /4)(fcf-A+ 2 /4) 

. A-l 

/ dxi dx 2 d5id5 2 d 2 x± Y\ d 3 Xit 



I 2 j j ax 1 ax 2 dd 1 dd 2 crx± |j_ a 
x $*({xr +6i/2,x x },{Xi +8 2 /2,x ± },...)<P({xi -6^/2,x ± },{x2 - 5 2 /2, x^}, . . .) 

1 f dy~ e iy~{x A p++l + ) e iu-{x B p + -l + ) 

x - / du~ 



2 J 2tt (k+ - ie)(l+ - ie) 

x (h- A/2\^(0,x x ) 1 +'iP(y-,x x )\k 1 +A/2)(k 2 + A/2\F+(-u-/2,x x )F' 7+ (u-/2,x 1 _)\k 2 - A/2). (50) 

In the above equation, A + = k + + 1 + — xcp + , we continue to use A + to save writing. The only unphysical variables 
remaining are the two momenta k + ,l + . Again, following Ref. [3l|, we note that the the OFPD's also demonstrate 
minimal variation with A + and as a result that dependence may be ignored. Restricting attention solely on the part 
that depends on k + , l + we obtain the integrals 

dkdl e ik+ ( x 2 ~ x i) e il+ ^2 ~ x i ) e a+ _ "~ ) 
(2ni) 2 {k+ - ie){l+ - ie) ' 

This yields the obvious condition 9{x 2 — x^)0(x 2 — X\ + y~ — u~). Within the picture adopted by the convolution 
model, X\,x 2 are the locations of the centers of the two nucleons that contain the struck quark and the soft gluon 
which scatters off the struck quark, where as u, y are locations within a given nucleon. As a result in the case where 
the struck quark and the soft gluon originate within two separate nucleons (as is the case in this calculation), the 
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second 8- function is identical to the first and the overall integration over k + ,l + simply ycilds the physical condition 
that the nucleon containing the struck quark be situated spatially ahead of the nucleon containing the soft gluon. The 
contour integration over the complex space of l + sets their values to zero, as a result A + = xcP + ■ By inspection 
of Eq. (|28p . one notes that the only values that may be assumed by xc are or x^. As a result, we make the further 
approximation that xcp + << kf,k 2 and neglect xc in the denominators of the top line of Eg. (|50|) . Corrections 
to this approximation are either vanishing or suppressed by powers of x\. Following Ref. (29j |. and given that the 
nuclear wavefunction is peaked for values of ki+ — p + we make the approximation of replacing the off- forward parton 
distributions with the regular distributions, 

74 J ~ A/2|VX0,^) 7 + V >( Z y-,^)|fc 1 + A/2) = 2Q (^^, ^ 2/ g (x) (51) 

r r U,- p iu~x B p + / ' Tn + A+ /I/ 2 _ 1,2 \ 

J ~2^ k+ {k2 + A / 2 l F - + (-^/ 2 ' ^)^ + (^/2, x ± )\k 2 A/2) = 2G (Jp - — , -jjt^) * V 9 ( x ) (52) 

It should be pointed out that the above approximations for the parton distribution functions have only been demon- 
strated to hold in convolution with nuclear density distributions which are peaked around the mean values of the 
momenta of the nucleons (see Ref. [29| for further details). 

Incorporation of the above simplifications leads to the following expression for the representative matrix element, 



T ~ / a 7 A 2 / dx-dx^dS^cPx^ H d 



3 r .piki&x „ik 2 S2 ix c p + (x 2 —x 1 ) 



i=3 

x $*({xf + x ± }, {xz + S£/2, n}, . . .)$({xi - 6i/2, x x }, {x 2 ~ 5^/2, x x }, ■ ■ ■) 

x 2irf q {x A )f g (x B ). (53) 

In the above equation, P ma x is the maximum allowed momentum in the nucleus, in principle this should be weakly 
dependent on the total number of nucleons in the nucleus, we will ignore this dependence in this effort. The integrals 

t- 1 



2vr 

are sharply peaked around the point 5 n = with a height proportional to P ma x ■ We thus replace these integrals with a 
saddle point approximation around the point S n = 0. All factors of the width of the integrals will be included within an 
overall unknown normalization constant B. This approximation simplifies the product of the nuclear wave-function 
and its complex conjugate &*({x^[ + 5^/2, x±}, {x^ + S 2 /2, x±}, . . .)<f>({x^ — 8^/2, x±}, {x^ — ^/2, x±}, ...)—* 
^*({ x i i x ±}i i x 2 > x -i-}i ■ ■ -)^{{ x i > x -l}i { x 2 i x ±-}i ■ ■ ■) — p[( x ii x ±)i ( x 2, x ±)]- This is the two nucleon correlation func- 
tion in a nucleus. We make the last approximation of replacing this with a product of two single nucleon densities 
and an unknown constant D. i.e., 



p(x 1 ,x 2 ) = Dp{x l )p{x 2 ). (54) 

this leads to the final simplified expression for the representative matrix elements which has already appeared previ- 
ously in Refs. 0, HU , 



T = j dxidxz d 2 x ± Cp(x]; ,x ± )p(x^ ,x ± )e ixGP+ ^- x ^6(x^ - x^)A 2 f q (x A )x B f g {x B ), (55) 

where all normalization constants including factors of 2tt and p + have been absorbed into the dimensionful constant C . 

Such a constant may not be determined from first principles in our approach, but should, in principle, depend on 
the kinematics of a given experiment. As our derivation has demonstrated, it also may have a mild variation on the 
nucleus chosen. However, we will assume that the overall constant be independent of the nucleus chosen. Ostensibly, 
the overall constant also has no dependence on the fragmentation of the escaping jet and thus may not depend on 
the momentum fractions of the detected hadrons Zi,z 2 nor on the number of hadrons detected. Once a data point 
for a choice of momentum fraction (fractions) of a detected hadron (hadrons) is described by tuning C, it is set for 



16 



all other momentum fractions and nuclei that may be used in the given experimental kinematics. The constant has a 
strong dependence on p + , however, as in most experimental setups the energy per nucleon in an accelerated nucleus 
is held fixed for different nuclei, there is no variation across nuclear targets for different total momenta . 

In the present section, we have introduced a number of approximations in the evaluation of the medium modification 
factor outlined in Eq. (|28|) . The approximations have resulted in a very much simplfied expression for the modification 
factor. In the subsequent section we will evaluate the modification of the dihadron fragmentation functions using the 
expression from Eq. (|55[) . Such a computation will require a certain model of the nuclear density distribution p{x). 
Calculations for two separate distributions (gaussian, hard-sphere) will be presented. 



IV. NUMERICAL RESULTS 



The modification of single and double hadron fragmentation functions in nuclei has been studied experimentally by 
the HERMES experiment at DESY 0, [2l|, [HI ■ The modification is presented as a ratio of the fragmentation of a hard 
parton produced in the DIS off a large nucleus, versus a deuterium nucleus. Within the kinematics of the experiment, 
the DIS off the deuterium nucleus leads to minimal modification of the fragmentation process. Fragmentation of the 
jet produced in such a process is assumed to be identical to fragmentation in vacuum. 

In the rest frame of the struck nucleus the experiment can measure both the forward momenta of the produced 
particles Pi ,p 2 , the forward momentum of the virtual photon q~ = v as well as the virtuality of the photon Q 2 . The 
momentum fractions are thus given simply as the ratios Z\ = V\ l v an d z 2 = V 2 l v - The modification of the single 
fragmentation function at a momentum fraction z is represented by the ratio of the multiplicity of hadrons N(z)dz 
with a momentum fraction in the range from z to z + dz produced in DIS off a large nucleus versus that produced in 
DIS off a deuterium nucleus: 

Theoretically this ratio is equated with the ratio of the medium modified fragmentation function versus that in the 
vacuum, i.e., 

D(z,v,Q 2 ,A) 

Rl{z) = D(z,Q 2 ) ■ (5?) 

Mcsurements for a Nitrogen (N) and Krypton (Kr) nucleus are presented as the square and circular points in Figs. 16181 
The medium modification of the associated hadron fragmentation function is obtained by measuring the number 

of events with at least two hadrons with momentum fractions z\, z 2 (N\ (zi, Z2)), and the number of events with at 
least one hadron with momentum fraction z (Na(z)). As in the case of the single fragmentation function a double 
ratio is presented using similar measurements off deuterium, 



y 

N A (z) 

R 2 {z 2 ) = -. (58) 

This double ratio for different z 2 is plotted in Figs. [7] and [9] as the square points for nitrogen and the circular points 
for Krypton. In the limit of low multiplicity per event, and the assumption that there is minimal modification in 
Deuterium the above ratio is theoretically estimated as 

f 1- * 2 dz 1 D„(z 1 ,z 2 ,u,Q 2 ,A) 
J Q-5 : 

f dzD„(z,u,Q 2 ,A) 

MZ2) = J™, — • (59) 

J os ^dz^jz^Q 2 ) 
£ 5 dzD q (z,Q2) 

In the remainder of this section we focus on the evaluation of the medium modification factor and the medium 
modification of single and double fragmentation functions. In the proceeding section, a set of approximations to the 
medium modification factor of Eq. (f!?5|) were carried out. There exist four terms in Eq. (f25)) . each of which may be 
expressed in general as in Eq. (TM|) . The generic term, under the approximations instituted in the previous section 
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may be expressed as Eq. (f55j) . Substituting Eq. 
Eq. (j2"5)) leads to the expression, 



55|) back into the complete expression for the modification factor 



T^ (x,xl) = / dx 1 dx 2 d 2 x±Cp(x 1 , x±)p(x 2 , dx±)8(x 2 — x 1 )A 2 



f q (x + x L )x T f g (x T )^l-e-™-P + ^~ x i^+f q (x)x L f g (x L )^ 



(60) 



which contains all the four terms of Eq. (|28|) . To proceed further, an expression for the nuclear density distribution 
has to be substituted in the above equation. Two different choices of such a distribution will be studied in this 
article: a Gaussian distribution, which essentially concentrates most of the nucleons towards the center of the nucleus 
and a hard sphere distribution where the nucleons are distributed evenly over the entire nuclear volume. No doubt, 
most nuclei lie somewhere in between these two extremes. We refrain from using a Woods-Saxon distribution [32| . 
as it becomes analytically intractable. The hard-sphere distribution is a good approximation of the Woods-Saxon 
distribution for large nuclei. In the following we present results for each of the different cases in turn. 

The computation of the modification factor will be followed by the calculation of the modification of the frag- 
mentation functions in Eq. (|26[) for the dihadron fragmentation function and for the single fragmentation function 
following Ref. jij. As is clear from Eq. (|26|) . the evaluation of the medium modified fragmentation function also 
requires the input of a vacuum dihadron fragmentation function and vacuum single fragmentation functions. The 
vacuum dihadron fragmentation functions may not be predicted entirely with QCD but require a measurement at a 
given energy scale. In the absence of experimental data we follow Ref. [l2[ and parametrize two-particle correlations 
in the tuned Monte-Carlo event generator JETSET The vacuum single fragmentation functions are tabulated 
and parametrized by various collaborations. We choose the parametrization of Ref. [3~4| . This parametrization is also 
consistent with single particle distribution from our chosen event generator. 



A. Gaussian distribution 



Imagine a large nucleus with A nucleons at rest, at a given time t . The mass of each of the nucleons is M, and 
the nuclear radius is Ra- The spatial distribution of nucleons is assumed to be that of a Gaussian, i.e., 



P(y,x±) = A) exp 



2R\ 





x\ - 




2 

y ] 


= Po exp 


- 2R\ . 


exp 


- 2R\ \ 



(61) 



where po is the normalization factor, y is the coordinate in the 3-direction and x± is the two dimensional vector 
orthogonal to it. The lightcone vector is y~ = to — y at to = 0. The density is normalised to unity, i.e., 



J dy-dWpoe"^ eH = 1 =* Po = (62) 

We then boost the nucleus in the y-direction to a frame where its momentum Pa = A[\J M 2 + p 2 , 0, 0,p] is very 
large. For a large enough boost, we may approximate the light cone coodinates in the boosted frame y^ = y~ /2j, 
where 7 = p + /AI. The intergral of the density in the coordinates of the boosted frame becomes, (dropping the 
subscript b), 



J dy-d 2 x JL 2-fp(y-, x 1 _) = J ' dy~ d 2 'x ± 2 7 p e e 2(«a/^> 2 = 1 (63) 

The overall factor of 27 may be absorbed into a redefinition of po. This expression for the density is then substituted 
into Eq. (|60[) . There exist two kinds of terms in Eq. (|60[) : those with and without the phase factor exp[±ix lP + {x 2 — 
Xi)]. Summing the terms without the phase factor, one obtains, 



T$ x \x,x L ) = J ' dxldx-d'x^C- 



'l" 



YR\ 



x-)A 2 



x p e R Ae 2 ( r a/^) 2 f g (x + x L )x T fg(x T ) 



(64) 
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It is now trivial to carry out the Gaussian integrals over the position variables which yields factors of 

2tt, 7 and R A . Each power of Ra yields a factor of A 1 / 3 . Counting the power of A and absorbing all other factors of 
p + and 2tt into the overall normalization constant, we obtain the form as in Ref. Q 



T£W(x,x L ) = CMR A ff(x + x L )x T f g (x T ), (65) 
where we have approximated the quark structure function in a nucleus to be simply A times that in a nucleon 



ff{x + x L ) = Af q (x + x L ). 



The terms with the phase factor undergo a similar intergration, these may be summed to obtain 

Tqg {1 \x,X L ) = - J dx^dx^Sxt 2C ^y R Q 0( x 2 ~^i)A 2 



x p e r ac ^a/^)- f q ( x + x L )x T fg{x T )cos[xLP [x 2 - x i )] ( 66 ) 

The part of the above exponentials that depend on the positions of the two nucleons may now be integrated out as 
before. Adding the the contributions (both with and without the phase factor) results in the final expression for the 
medium modification function, in the Gaussian density approximation, 



(l_ e -*i/**) 5 ( 67 ) 

where xa — ^/(MRa)- The above expression for the modification factor is then used in Eq. (|2~7| to obtain the modified 
splitting functions and the medium modified fragmentation functions in the Gaussian density approximation. There 
still remains the overall factor C ~ C[xfg (x)] which represents the correlation between the nucleons and the gluon 
distribution at small x. We will consider this factor a fit parameter. It is determined for one measurement of the 
fragmentation function (single or double) at one value of momentum fraction of the hadrons in the DIS off one nucleus. 
We will use that value to predict medium modification of both single and dihadron fragmentation functions. 

The modifification of the single and double fragmentation functions in the DIS off nuclei {D(z), D(z\, 22)) has been 
measured by the HERMES experiment at DESY. The modification for the single inclusive distribution is plotted 
in Fig. [5] for Nitrogen (squared points) and Krypton (circular points) targets . These are then divided by the 
identical measurement on the deuterium nucleus, assuming minimal modification of the fragmentation function in 
deuterium. One notes a continuing suppression with increasing momentum fraction z and increasing nuclear size. The 
suppression factor is expressed in our calculation as the ratio of the medium modified to the vacuum fragmentation 
functions. The value of C is set by fitting the point at the lowest z for Nitrogen targets. Three different choices for 
C = 0.006, 0.007, 0.008 are presented. Arguably the best fit is with C = 0.006. The variation with increasing z is 
then a prediction from Eq. (|26|) and is shown as the solid line for the N nucleus. The calculation of the suppression 
factor in Kr requires nothing more than a change of the number A from 14 to 81. One notes good agreement with 
the data. However the agreement seems to deteriorate for smaller z. 

The modification of the dihadron fragmentation function as a function of momentum fraction for TV is presented 
in Fig. [51 It should be stressed once again that there are no free parameters in the calculation of the dihadron 
fragmentation function. The sole parameter C was set by comparison to the medium modified single fragmentation 
function (Fig. [5]). It turns out that within the choice of parameters the comparison with the TV nucleus is very good. 
The modified dihadron fragmentation function for Kr is shown in Fig. ((7]). One notes the agreement deteriorates for 
increasing z. The best fit to the data is obtained, once again with C = 0.006. 

It would seem that the Gaussian approximation for the nuclear density distribution provides a good fit to the data 
for both single and double pion distributions in the case of N. The comparison is not so good in the case of Kr. 
However, it should be noted, that the Gaussian density approximation, which is close to the actual density distribution 
in N is far from the actual nuclear density distribution in a larger nucleus such as Kr. In the next subsection, we 
will use a hard sphere distribution which is a better approximation for large nuclei. 



C_ 

X A 



ff(x)x T f?(xr) + ff(x)(x L + x T )f^(x L + X T ) 



B. Hard-sphere distribution 



In this subsection nuclear density distribution will be assumed to be a hard-sphere, 

p(y, x±) = p o 0{Ra - Jy 2 + x\). 



(68) 
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FIG. 6: The medium modification of the single fragmentation functions compared to data from the HERMES collaboration. A 
Gaussian distribution of nucleons in a nucleus is used. The dotted lines indicate a C = 0.006, the dot-dashed lines a, C — 0.007 
and the solid line a C = 0.008. 




FIG. 7: Results of the medium modification of the quark fragmentation function in cold nuclear medium in a N and Kr 
nucleus versus the momentum fraction of the associated hadron. The momentum fraction of the trigger, is held above 0.5 and 
integrated over all allowed values. A Gaussian distribution of nucleons in a nucleus is used. See text for details. 

The normalisation, 

J dyd 2 x±p(y, x±) = po^-R\ = 1, (69) 

implies pq = 3/(4irR\). As before, we now boost the nucleus to a frame where its forward momentum fraction is 
Ap + . The nuclear density in terms of the variables y~ , x± is given as 

p(y~, x±) = 2 1Po e(R A - sJy-^Arf + x\). (70) 
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We undertake the approximation that (xl +%T)f£ f (%L +%t) — xtJ^{xt) for the gluon densities in Eq. (f60|) . This 
leads to the simplified expression, 

T*(x,x L ) = f dx^dx 2 d 2 x 1 _2C 18 2 ? ' flfe - x^)A 2 Po 9(R A - \JxfA^ + x\) 



x p Q 9(R A -yx 2 47 2 +x 2 L )f q (x + x L )x T f g (x T )[l -cos{x L p + (x 2 - x x )}]. (71) 
Integrating over the transverse coordinate, one obtains 

rRA/2-y rRA/2~< 2 

Tq g {x,x L ) = dx^ dx 2 \R\ - (^2 7 ) ] 

J-R A /2-t Jx~ 

x CA^/,f (x + x L ).T T / 9 (a;T)[l - cos^p+fe - x^)}]. (72) 

In the above equation, one of the factors of A has been absorbed into the definition of the quark distribution function 
in a nucleus fa( x ) = Af q (x). As in the case for the Gaussian density distribution, one may also absorb the gluon 
density XTfg(xr) into the over all normalization constant C ~ Cxrlgixr)- Integrating over X2 and x\ and absorbing 
all constants except for factors of mass number A, one obtains the final expresssion for the medium modification 
factor in the case of a hard-sphere distribution as 



T Q A Jx,x L ) = C 



X A 



1 X"a I XT 

- + 4-£ sin [ — 

X A 



*4< 



XL 
X A 



(73) 



Substituting the above result into the modified splitting functions in Eq. (|27[) . one can calculate the medium 
modified fragmentation functions for a hard-sphere nuclear density distribution. There is the overall constant C as in 
the case of the Gaussian density distribution. This is set by fitting to the data. Results for the medium modification of 
the single hadron fragmentation function in the hard sphere approximation are presented in Fig. [S] for three different 
values of the parameter C = 0.022, 0.024, 0.030 as compared with the identical data set as in Fig. [6] The data seem 
to prefer the value for C ~ 0.022. This fits both the lowest momentum fraction point in ./V and provides a somewhat 
better description of the large z region in Kr. In this sense a hard-sphere approximation seems to work better for the 
case of a large nucleus such as Kr. 




FIG. 8: The medium modification of the single fragmentation functions compared to data from the HERMES collaboration. 
A hard-sphere distribution of nucleons in a nucleus is used. The dotted lines indicate a C = 0.022, the dot-dashed lines a 
C = 0.026 and the solid line a C = 0.030. 
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The modification of the dihadron fragmentation function as a function of the momentum fraction in the hard- 
sphere density approximation for N is presented in Fig. Once again, there are absolutely no free parameters in 
the calculation of the dihadron fragmentation function. The sole parameter C was set by comparison to the medium 
modified single fragmentation function (Fig. [5]). It turns out that within the choice of parameters the comparison 
with the N nucleus is not as good as the case for the Gaussian density distribution. This is to be expected as the 
Gaussian distribution is indeed closer to the actual density distribution in N . 

The modified dihadron fragmentation function for Kr is presented in Fig. One notes that while the agreement 
with the data still deteriorates for increasing z, there exits an overall quantitative improvement in the fit over that 
obtained from the Gaussian density distribution. The best fit to the data is obtained, once again with C = 0.03. It 
should be pointed out, in passing, that besides the slight improvement in the results for the case of the Kr nucleus 
there remains no qualitative difference in the modification of the single and double fragmentation functions between 
the cases of a Gaussian distribution of nucleons and a hard sphere distribution. 




FIG. 9: Results of the medium modification of the quark fragmentation function in cold nuclear medium in a N and Kr 
nucleus versus the momentum fraction of the associated hadron. The momentum fraction of the trigger, is held above 0.5 and 
integrated over all allowed values. A hard-sphere distribution of nucleons in a nucleus is used. See text for details. 



The curves presented in Figs. are the main results of this paper. They demonstrate that including next-to- 
lcading twist corrections in the DIS off nuclei allow for a better understanding of the measured data. However it is 
clear that one needs to go further to understand the behaviour displayed by the one and two particle distributions 
in the DIS off large nuclei such as Kr. For smaller nuclei such as N the next-to-leading twist contributions seem to 
provide a satisfactory account of the experimental behaviour. 



V. DISCUSSIONS AND CONCLUSIONS 



The focus of this paper has been on the medium modification of dihadron fragmentation functions in the semi- 
inclusive DIS off a large 1 nucleus. We have first generalized the formalism for the modification of the single fragmenta- 
tion function in a dense medium to the modification of the dihadron fragmentation function. The modification arises 
from the inclusion of next-to-leading twist contributions to the double differential inclusive cross section for observing 
two hadrons within a jet produced via leptoproduction from a nucleus. Higher twist contributions are suppressed 
by powers of the hard scale Q 2 and are thus ignored in the DIS off a nucleon target. A class of these higher twist 
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contributions are however enhanced by nuclear size A 1 ! 3 (23| and lead to the medium modification of both single and 
dihadron fragmentation functions. 

We have demonstrated (in Sec. II) that the medium modification of the dihadron fragmentation functions may 
be simply computed as a sum of convolutions involving medium modified splitting functions and vacuum dihadron 
fragmentation functions and and a qualitatively new contribution involving a new modified splitting function P q ^ qg 
and a product of quark and gluon single hadron fragmentation functions [Eq. I|26j) ]. Two of the modified splitting 
functions were shown to be identical to that in the case of the single fragmentation functions (i.e., AP q —> qg and 
AP q ^ gq ). The new splitting function AP q ^ qg was shown to be equivalent to AP q —> qg but without any virtual 
correction. Thus all modified splitting functions depend on the same medium modification factor T gg [Eq. I28| . 
identical to the case of the single fragmentation functions. 

These higher-twist contributions to the semi-inclusive DIS cross section involves multiple scattering and gluon 
bremsstrahlung encountered by the struck quark and includes the well known Landau-Pomeranchuck-Migdal (LPM) 
interference effect. This interference effect leads to the appearance of the off-forward parton distribution functions 
(OFPD) [3(| in the nuclear modification factor T^ g . Based on phenomenological models of OFPD's we approximated 
these by the regular forward parton distributions. The final result depends on the density distribution of the nucleons 
in the nucleus [Eq. ([60]) ]. Results for two different distributions were presented: Gaussian (Figs. Ed) and hard-sphere 
(Figs. I8l9p . The two distributions demonstrated qualitatively similar effects on both single and double inclusive 
spectra. 

The evaluation of the modification factor involved an overall normalization constant which represents the correlation 
between the nucleon containing the struck quark and the nucleon containing the soft gluon which scatter off the struck 
quark. This was set by fitting to the overall experimental data on single inclusive distribution for DIS off N . The 
computation of the modification of the single fragmentation function on Kr and the dihadron fragmentation function 
on both N and Kr thus involves no free parameters. The comparison of the theoretical prediction for the medium 
modified fragmentation functions (both single and double) with the experimental data for DIS off a N nucleus is very 
good. The comparion for the case of the Kr nucleus is somewhat satisfactory. A possible cause for this may lie in the 
inclusion of only the twist four contributions. These are contributions proportional to A 1 / 3 /Q 2 that are suppressed 
by a power of Q 2 yet enhanced by a factor of A 1 / 3 . Physically this means that the struck quark may undergo at 
most two more scatterings off soft gluons prior to exiting the nucleus and fragmenting. Further scatterings necessarily 
involve higher twist parton correlation functions and require the inclusion of further powers of A 1 / 3 /Q 2 . While, such 
contributions may be unimportant in the case of N, they may become necessary for the case of Kr and heavier nuclei. 
A systematic inclusion of all orders of A 1 / 3 /Q 2 in the modification of the fragmentation process will involve a far 
more complicated calculation and is beyond the scope of this paper. 
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